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Abstract. We examine the energy surfaces of the driven Rabi model, also known 
as the biased or generalised Rabi model, as a function of the coupling strength 
and the driving term. The energy surfaces are plotted numerically from the known 
analytic solution. The resulting energy landscape consists of an infinite stack of sheets 
connected by conical intersection points located at the degenerate Juddian points in the 
eigenspectrum. Trajectories encircling these points are expected to exhibit a nonzero 
geometric phase. 


The driven, biased or generalised Rabi model has hamiltonian PEIE] 

H = oj a^a + g ax{a^ + a) + A az + e ax , (1) 

where ax and az are Pauli matrices for a two-level system with level splitting A. The 
interaction between the spin and the single-mode bosonic held of frequency oj is via the 
coupling g. The bosonic creation and destruction operators and a satisfy [a, = 1. 

The Rabi model is well known as arguably the simplest model for light interacting with 
matter. As such, it has a long history in quantum optics||] The addition of the driving 
or bias term eax breaks a Z 2 symmetry (parity). This additional term allows tunnelling 
between the two atomic states. The driven Rabi model ([^ is relevant to the description 
of various hybrid mechanical systems W- In particular, of a micromechanical resonator 
via coupling to a Cooper-pair box [SI ElIZl El- 

Using the analytic solution obtained by Braak [1] for the energy eigenspectrum, 
we explore the energy levels of this model as a function of the parameters g and e. 
Specihcally, the Nth eigenvalue is given by Ej^ = — g‘^/oj, where is the Nth zero 

of 


G,{x) = A‘^R+{x)R-{x) - R^{x)R-{x ), 


( 2 ) 


f For realisations in ion traps and in both cavity and circuit quantum electrodynamics, the reader is 
referred to Ref. |1] and references therein. 
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Figure 1. Energy eigenspectrum of the driven Rabi model as a function of the coupling 
g for parameter values A = 0.7 and a; = 1. From left to right, e = 0, ^ and 1. 


where 


n=0 



n=0 


X — nov ± e 



n 


( 3 ) 

( 4 ) 


K^{x) is defined recursively by nKf; = — K^_2 with initial conditions 

Kq = = fo{x), and 



^ + 1 

u 2g 


— X ± e + 




X — no; ± e 


( 5 ) 


The function G±{x) is obtained as a consistency condition for two different series 
expansions in a representation of the bosonic creation and annihilation operators in the 
Bargmann space of analytical functions [1]. It can also be derived using the extended 
coherent states approach |2] and written in terms of confluent Heun functions. The 
eigenstates have also been obtained in terms of confluent Heun functions PE]. 

The lowest energy levels are shown in Figure 1 as a function of g for fixed A and 
e. For e = 0 this is the well known plot featuring the Juddian crossing points (TU] for 
which Kus [11] provided a proof that for each value of N there are N crossings if A is in 
the range 0 < A/o; < 1. More generally for k < A/uj < k + 1 there are N — k crossing 
points m- There are also such crossing points in the eigenspectrum when e is an integer 
multiple of HP- These are shown in Figure 1 for e = | and e = 1. It has also been 
argued that for a given value of N there are N level crossings for 0 < A/w < \/l + 2e/uj, 
reducing to N — k crossing points for A in the range jl20] 

\/k‘^ + 2ke/oj < A/oj < \/(/c + 1)^ + 2{k + l)e/a;. (6) 


§ The eigenspectrum is symmetric in e, so this relation holds in general with e replaced by |e|. 
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Figure 2. Energy landscape of the driven Rabi model as a function of the parameters 
g and e for the values A = 0.7 and w = 1. The landscape consists of a series of discrete 
sheets labelled by integer N and connected by conical intersection points for A > 1. 


Here we explore the energy landscape as a function of the parameters g and e for 
hxed A. As we shall see, the effect of varying the parameter e is to induce conical 
intersection points. The tips of the cones are precisely the degenerate Juddian points. 
Figure 2 shows the structure of the interconnected energy surfaces as a function of g 
and e for A = 0.7. This hgure should be viewed in conjunction with Figure 1, which 
gives the cross-sections at e = 0, ±^,±1 of the energy landscape shown in Figure 2. 
This figure also shows crossing points for = 0 as the parameter e is varied. These are 
determined by the energies 

E = Nu± VA2 + e2 . (7) 

Two more rehned views of the energy surfaces are shown in Figure 3jj[| It is clear that the 
energy surfaces for g > are connected by isolated conical intersection points. These 
are shown further in Figure 4, which illustrates the elementary cones. 

The overall structure of the energy landscape is an inhnite stack of sheets labelled 
by integers N = 0,1,2,... connected for g > Ohy conical intersection points for A > 1. 

II Note that the overall structure of these figures can be made simpler by plotting E + g'^foj on the 
z-axis which has the effect of turning the paraboloids into planes. 
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Figure 3. Close up views of the energy landscape depicted in Figure 2 showing the 
energy surfaces connected by conical intersection points. 




Figure 4. Elementary cones in the energy spectrum for A = 0.7 with uj = 1. The plot 
on the left shows the “lowest” cones at e = 0, ^ and 1. Also shown is a magnification 
of the cone centred at e = 0. 
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The conical intersection points are the exceptional points in the eigenspectrnm which 
are located on lines in the planes defined by e = ±^nu, n = 0,1, 2,.... These lines are 
the intersections of the snrfaces E = Niu — g‘^/uj — e and E = N 2 U — g‘^/u + e for different 
^ 1 , ^2 > 1 with A^i — A ^2 = ^e/co. The precise locations of the conical intersection points 
along these lines are determined by the solntions of the constraint relations [SI IT2] . 
Their number is dependent on the value of A/u. From the above consideration of the 
number of crossing points, the number of conical intersection points can be tabulated 
for given A/u and e/u. In this sense the precise geometry of the energy landscape is 
thus known. For example, for the value A = 0.7 considered in the figures, there are N 
conical intersection points for a given value of N located in each of the planes where 
e = A^nu, n = 0,1, 2,.... Figure 4 identifies the most elementary of the cones, with a 
magnification of the cone centred at e = 0. From the general result ([^ for the number 
of crossing points we know that this particular cone vanishes for A/uj > 1. 

Trajectories encircling conical intersection points are expected to exhibit a nonzero 
geometric phase. Geometric phases uni HU have been investigated in the context of 
the Rabi model for some time, and not without debate (see, e.g.. Refs. Ha HHl HU 
and references therein). For systems like the Rabi model the Berry phase is induced by 
integration over a unitary transformation parameter ip with the angle variable ip varying 
slowly from 0 to 2n [TH]. Here we have demonstrated conical intersection points in the 
energy landscape for the driven Rabi model by varying the system parameters g and e. 
Since analytic expressions for the corresponding eigenstates are known for this model, 
the geometric phases obtained by integrating trajectories encircling the conical points 
should, in principle, be calculable. 

As a final remark, we note that the original justification for studying Juddian 
points was that, because they were exact solutions, they could be useful in testing 
and improving various approximation schemes. Here, using the analytic solution of the 
driven Rabi model, we have demonstrated that the degenerate Juddian points are the 
conical intersection points in the energy landscape. Degenerate Juddian points exist in 
a range of models related to the Rabi model. It is hoped that these models can be added 
to the different contexts in which conical intersection points have been observedjj] The 
conical intersection points are expected to play a crucial role in the dynamics of the 
driven Rabi model, which is as yet only partially [3] explored. 
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